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A stochastic gravitational wave background causes the apparent positions of distant sources to 
fluctuate, with angular deflections of order the characteristic strain amplitude of the gravitational 
waves. These fluctuations may be detectable with high precision astrometry, as first suggested by 
' Braginsky et al. in 1990. Several researchers have made order of magnitude estimates of the upper 

| limits obtainable on the gravitational wave spectrum f2 gw (/), at frequencies of order / ~ lyr -1 , 

both for the future space-based optical interferometry missions GAIA and SIM, and for VLBI 
interferometry in radio wavelengths with the SKA. For GAIA, tracking TV ~ 10 6 quasars over 
a time of T ~ 1 yr with an angular accuracy of A9 ~ lO^as would yield a sensitivity level of 
f2 gw ~ (A9) 2 / (NT 2 H 2 ) ~ 10 -6 , which would be comparable with pulsar timing. 
• In this paper we take a first step toward firming up these estimates by computing in detail the 

statistical properties of the angular deflections caused by a stochastic background. We compute 
analytically the two point correlation function of the deflections on the sphere, and the spectrum as 
a function of frequency and angular scale. The fluctuations are concentrated at low frequencies (for 
f ) ' a scale invariant stochastic background), and at large angular scales, starting with the quadrupole. 

The magnetic-type and electric-type pieces of the fluctuations have equal amounts of power. 
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, A. The Stochastic Gravitational Wave Background 

There is great interest in detecting or constraining the strength of stochastic gravitational waves (GWs) that may 
have been produced by a variety of processes in the early Universe, including inflation. The strength of the waves is pa- 
J^j . rameterized by their energy density per unit logarithmic frequency divided by the critical energy density, f2 gw (/). Cur- 
rent observational upper limits include (i) the constraint fi gw < 10- 13 (//10- 16 Hz)" 2 for 1CT 17 Hz < / < Kr 16 Hz 
, from large angular scale fluctuations in the cosmic microwave background temperature [l[ ; (ii) the cosmological nucle- 
■ osynthesis and cosmic microwave background constraint J din / f2 gw (/) < 10~ 5 , where the integral is over frequencies 
rj< ■ / ~ 10_15 Hz 0; ( Hi ) the P ulsar timing limit fl gw < 10" 8 at 1CT 9 Hz < / < 1(T 8 Hz 0; (iv) the current LIGO/VIRGO 
upper limit ft gw < 7 x lfr 6 at / - 100 Hz Q; and (v) the limit J din f fi gw < 10" 1 for 10" 17 Hz < / < 10~ 9 Hz from 
VLBI radio astrometry of quasars. 
7-H . Many new techniques also promise future measurements of these primordial GWs. Firstly, it has been shown that 
such a GW background would leave a detectable signature in the polarization of the cosmic microwave background 
[CMB;[^,Q, which will be measured by many current and future observational efforts 0-Q11- The planned space-based 
interferometer LISA will also set limits on the primordial stochastic gravitational wave background (SGWB) [ijj]. The 
planned successor to LI SA, the Big Bang Observer, is a space- based interferometer mission designed primarily to detect 
the primordial SGWB 18]. Finally, Seto and Cooray have suggested that measurements of the anisotropy of time 
variations of redshifts of distant sources could provide constraints of order f2 gw < 10 -5 at / ~ 10~ 12 Hz [19]. For 
more details on GWs, the search for them, and the SGWB, see the review articles (ll I20I l2l| . 



INTRODUCTION AND SUMMARY OF RESULTS 



B. High Precision Astrometry 

The possibility of using high precision astrometry to detect GWs has been considered by many authors [221433 ]. 
There was an early suggestion by Fakir (2(| that GW bursts from localized sources could be detectable by the angular 
deflection A9 to light rays that they would produce. Fakir claimed that A9 oc 1/b, where b is the impact parameter. 
This claim was shown later to be false, and in fact the deflection scales as 1/b 3 [2ol W\ . Therefore the prospects for 
using astrometry to detect waves from localized sources are not promising [341 ] . 

However, the situation is different for a SGWB, as first discussed by Braginsky et al. [24j. For a light ray propagating 
through a SGWB, one might expect the direction of the ray to undergo a random walk, with the deflection angle 



growing as the square root of distance. However, this is not the case; the deflection angle is always of order the 
strain amplitude h rms of the GWs, and does not grow with distance [i] 0,0, HH- Specifically, a SGWB will cause 
apparent angular deflections which are correlated over the sky and which vary randomly with time, with an rms 
deflection <5 rms (/) per unit logarithmic frequency interval of (see Eq. [3] below) 



<W/) ~ W/) ~ y^gw(Z). (1) 

Suppose now that we monitor the position of N sources in the sky, with an angular accuracy of AO, over a time 
T. For a single source, one could detect an angular velocity (proper motion) of order ~ A8/T, and for N sources, 
a correlated angular velocity of order ~ A0/(T\/~N) should be detectable. The rms angular velocity from (Q~]) is 
w rms (/) ~ / S rms (f) ~ Ho ^/Ogw(/)> and it follows that one should obtain an upper limit on fi gw of order [33| 

AO 2 

<W/) < ^ (2) 

This bound will apply at a frequency of order / ~ 1/T. It will also apply at lower frequencies [33| since the 
angular velocity fluctuations are white (equal contributions from all frequency scales) , assuming a flat GW spectrum 
f2 gw = const. The quantity that will be constrained by observations is roughly this total fi gw , J'- <T _ 1 d\n ffl gw (f). 

The advent of microarcsecond astrometry has started to make the prospects for constraining GW backgrounds more 
interesting. The future astrometry mission GAIA (Global Astrometric Interferometer for Astrophysics) is expected to 
measure positions, parallaxes and annual proper motions to better than 20 /ias for more than 50 million stars brighter 
than V ~ 16 mag and 500,000 quasars brighter than V ~ 20 mag [23]. Similarly the Space Interferometry Mission 
(SIM) is expected to achieve angular accuracies of order 10 uas. Estimates of the sensitivities of these missions to a 
SGWB, at the O gw - 10" 3 - 10" 6 level, are given in Refs. [H HI OS . 

VLBI radio interferometry is another method that can be used to detect the astrometric effects of a SGWB on 
distant sources. This method detects the same pattern as that discussed in this paper for visible astrometry, and 
differs from astrometry using the GAIA satellite in its longer duration (tens of years versus a few years for GAIA), and 
in the smaller number of sources, on the order of hundreds, that have currently been measured using this method. In 
the radio, the planned Square Kilometer Array (SKA) is also expected to be able to localize sources to within ~ lOuas 
(35| . Jaffe has estimated that with 10 6 QSO sources, the SKA could achieve a sensitivity of order f2 gw ~ 10~ 6 |2a |. 

The astrometric signals due to a SGWB expected for a single object are quite small, on the order of 0.1 n&s yr , 
much smaller than the typical intrinsic proper motion of a star in our galaxy. We therefore propose to use quasars as 
our sources, since their extragalactic distances cause their expected intrinsic proper motions to be smaller than those 
expected from a SGWB [23|. The construction of a non-rotating reference frame using quasars in astrometric studies 
will remove the I = 1 dipolc component of the measured quasar proper motions, but will leave intact the I = 2 and 
higher multipoles which are expected to be excited by GWs. 

Using the estimate N ~ 10 6 (GAIA), AO ~ 10 /ias, T ~ 1 yr gives from Eq. © the estimate 



O gw < 10~ 6 

at / < 10 -8 Hz for astrometry. This is an interesting sensitivity level, roughly comparable with that obtainable with 
pulsar timing Q. 

Astrometry has already been applied to obtain upper limits on f2 gw using a number of different observations. 
First, Gwinn et. al analyzed limits on quasar proper motions obtained from VLBI astrometry, and obtained the 
upper limit f2 gw < 10 _1 for 10~ 17 Hz < / < 10 _9 Hz (22|. This limit was recently updated by Titov, Lambert and 
Gontier (36|. Finally, Linder analyzed observed galaxy correlation functions to obtain the limit f2 gw < 10~ 3 for 
10~ 16 Hz < / < 10- 10 Hz H3. 

All of these analyses used a relativly simple model of the effect of graitational waves on proper motions. In this 
paper we give a detailed computation of the spectrum of angular fluctuations produced by a stochastic background, 
including the relative strengths of E- and B-type multipoles for each order I. In a subsequent paper we will follow up 
with a derivation of the optimal data analysis method and a computation of the tt gw sensitivity level, to confirm the 
existing crude estimates of the sensitivity of future astrometric missions such as GAIA. 



[i] It is sometimes claimed in the literature that the deflection angle depends only on the GWs near the source and observer. In fact, this 
is not true, as we discuss in Appendix A. A similar claim about the frequency shift that is the target of pulsar timing searches for GWs 
is also false in general. 



C. Summary of results 



For a source in the direction n, the effect of the GW background is to produce an apparent angular deflection 
<5n(n,i). We first find a general formula for the angular deflection of a photon, for an arbitrary GW signal hy, 
emitted by a source that can be at a cosmological distance. This deflection is derived in Sees. ITT1 and Hill below, and 
is given by [cf. Eq. (|55|) ] 



Sn* =-} n j hij(0) - n i n , n k h jk (P) ~ ~ n l n k ) n j - 

• -2 f <S dC f dC"h jkfi (C")+n l [ C& dC' f C d(" (h jkd (C") + hki.AC") - h 3l . k (C)) 
Jo Jo Jo Jo 

Here, n is the direction to the source, wo is the emitted frequency of the photon, f parameterizes the path of the 
photon t(£) = to + woC, £ l (C) = —C^on 1 , /iij(r,x) is treated as a function of £ through this parameterization of the 
photon path, ( s is the value of £ at the emission event of the photon at the source, and the spacetime metric is 



ds 2 = a(r) 2 {-dr 2 + [Sy + hij(T,x)] dx*dx j } . 

We then specialize to the limit in which the sources are many gravitational wavelengths away and to plane waves 
propagating in the direction p to obtain a simple formula, which generalizes a previous result of Pyne et al [33j . We 
find that the deflection, as a function of time r and direction on the sky n, is given by 



_>•: 1 • p-m "" ' " 2' 



Sti 1 (t, n) = — | - — hjk(0)njn k - ~hij(Q)rij, 



where p is the direction of propagation of the GW, and hij(0) is the GW field evaluated at the observer, hij(r, 0). 

The main result of this paper is a computation of the statistical properties of the angular deflection resulting from 
a SGWB, which is carried out in Sees. IIVI and fVl The apparent angular deflection caused by such a GW background 
is a stationary, zero-mean, Gaussian random process. We compute the fluctuations in Sn by making two different 
approximations: (i) The GW modes which contribute to the deflection have wavelengths A which are short compared 
to the horizon size cHq 1 today, (ii) The mode wavelengths A are short compared to the distances to the sources; this 
same approximation is made in pulsar timing searches for GWs [38| . Since our calculations are only vaild for GWs 
with wavelengths much smaller than the horizon, the contribution from waves with wavelengths comparable to the 
horizon scale will cause a small deviation from our results (on the order of a few percent for a white GW spectrum). 

The total power in angular fluctuations is then 

(Sn(n,t) 2 ) = Cs = ^J dlnf f^j n m (f). (3) 
Taking a time derivative gives the spectrum of fluctuations of angular velocity or proper motion: 

(Sh(n,t) 2 ) = J d\nfH 2 n gw (f), 
which gives a rms angular velocity w rms (/) of order 



_ Q \ 1/2 

Wrms(/) ~ i?oV^gw ~ 10 Vasyr 1 ' 
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This is the signal that we hope to detect. 



We now discuss how the angular fluctuations are distributed on different angular scales, or equivalently how the 
power is distributed in the spherical harmonic index I. The total angular fluctuations can be written as 



TABLE I: First 10 nonzero multipole coefficients af E as defined in Eq. (|87[) an d plotted in Fig. [T] 



I 


E E 


2 


0.833333 


3 


0.116667 


4 


0.03 


5 


0.0104762 


6 


0.00442177 


7 


0.00212585 


8 


0.00112434 


9 


0.000639731 


10 


0.000385675 


11 


0.000243696 



(5n{n lt f) = / dlnfJ2 K»,i(f) 2 + d(/) 2 ] • (4) 

J 1=2 

Here 9f ms i{f) 2 is the total electric-type power in angular fluctuations per unit logarithmic frequency in multipole 
sector and 0f ms ;(/) 2 is the corresponding magnetic-type power. These quantities can be written as 



8? ms Af) 2 = 8Ls9Qv(f)a? Q , (5) 

where Q = E or B. The various factors in this formula are as follows. The factors gE and gs are the fractions of the 
total power carried by E-modes and B-modes respectively, and satisfy qe + 9b = 1- Their values are gE = .9s = 1/2, 
implying that electric and magnetic type fluctuations have equal power. The function a(f) describes how the power 
is distributed in frequency, and is the same for all multipoles, both electric and magnetic. It is normalized so that 
J d In /cr(/) = 1, and is given explicitly by [cf. Eq. (j3|) above] 



rr( f\ — / 2 ^gw(/) 

aUj /dln/'/'- 2 ^ g w(/')' W 

Finally, the angular spectra af E and af B describe how the power is distributed in different multipoles, starting with 
the quadrupole at I = 2, and are normalized so that 



$>? Q = 1 (7) 

1=2 

for Q — E and Q — B. We show that af E = af B , and this spectrum is plotted in Fig. [T]and tabulated in table 
HI These coefficients are well fit by the power law af E — 32.34 l~ 4 - 921 . We note that the result for the quadrupole, 
a EE — 5/6, has previously been derived using a different method in Ref. |33j . 



II. CALCULATION OF ASTROMETRIC DEFLECTION IN A MINKOWSKI BACKGROUND 

SPACETIME 



A. Setting the stage— Minkowski Calculation 



We will first calculate the angular deflection due to a small GW perturbation on a flat background metric, 



ds = g^ydx 11 dx v — —dt + (Sij + hij)dx l dx 3 . 



(8) 




FIG. 1: Here we plot the coefficients af E as defined in Eq. (|87[) vs - multipole 



We are considering the effect of these GWs on a photon traveling from a source to an observer, with an unperturbed 
worldline Xq(X) = wo(A, — An) + (to, 0,0,0), where — n is the direction of the photon's travel, wo is its unperturbed 
frequency, and the photon is observed at the origin at time to- The photon's unperturbed 4- momentum is given by 
kg = wo(l,-n). 

To calculate the geodesies that the photon, source and observer follow, we need the connection coefficients in this 
metric. There are three non-zero connection coefficients: 



First, using the geodesic equation 



,00,-y 



(10) 



it is straightforward to verify that the paths of stationary observers in these coordinates are geodesies. Therefore we 
can assume that both the source and observer are stationary in these coordinates, with 



c* (t) = x\ = constant. 



The afhne parameter of the source is therefore 



UJQ 



B. Photon Geodesic 



Next, we solve the geodesic equation (flOl) for the path of a photon traveling from the source to the observer in the 
perturbed metric. We write this path as the sum of contributions of zeroth and first order in h, 

x a (X)=x%(X)+x<?(X). (11) 

Similarly, the photon 4-momentum is 

k a {\) = A£(A) + fc?(A). (12) 

We note that the connection coefficients are all first order in h, so keeping only first order terms, we will use only the 
unperturbed photon 4-momentum in the geodesic equation, yielding 

M O 2 

u x l i 71, /i o\ 

-^2 = — J n *■ ' 

-jTif = --^[-Irthufl + n l n j (h kiJ + h kjA - h ijjk )}. (14) 

We now integrate the geodesic equation to obtain the perturbed photon 4-momentum and trajectory. The right 
hand sides are evaluated along the photon's unperturbed path from A = at the present time back to A, since they 
are already first order in h. We define 

X ij (X) = f dX'h ijfi {X'), JvkW= I dX'h^kiX'), 

Jo Jo 

Kii{X) = [ dX' [ dX"h ij>0 (X"), C ijk (X) = [ dX' [ dX"h ij , k {X"), (15) 

Jo Jo Jo Jo 

where hij(X) means hij(to + i^oA, — woAn). We find 



4W = -^nV/C y (A) + IoX + K , 4(X) = -^-n l S tJ + J*X + L' Q , (16) 

where Iq, Jq, Kq and L J are constants of integration, and we have defined the quantities 

Rii (A) = [-21 t3 (A) + n k (J ijk (A) + J jki (A) - J ikj (A))] , (17) 



2 



Sij (A) = [-2Kii (A) + n k (C ijk (A) + C jki (A) - C ikj (A))] . (18) 
C. Boundary conditions 

We determine the eight integration constants I , Jq, K and L J using the boundary conditions of the problem, 
namely that the photon path passes through the detection event x^ bs — (to, 0,0,0), that it is null, that the photon 
is emitted with the unperturbed frequency uq, and that the photon path intersects the path of the source at some 
earlier time. 



1. Photon path must pass through detection event 

First, the perturbed photon trajectory must pass through the detection event t = to, x 1 = 0. Therefore, 



^(0)=^(Q)+<(0) = (t , 0,0,0), 

giving 



K = ^W/C y (0) = 0, L = y n'Si, (0) - 0, (19) 

where we have used the fact that by definition /Cy(0) = 5^ (0) = 0. 
2. Photon geodesic is null 

The geodesic of the photon must be null, which gives one more constraint: g^k^k^ — 0. This is already true 
to zeroth order. To first order we get: 

= h^k^k^ + 2rj llv k^kQ , 

where fcg = uo(l,— n). Inserting the expression for the perturbed 4-momentum k" given by Eqs. (|15[) , ([1 
and p7|) , and simplifying using 



_d_ 

dA 



hij — uj hij :0 — ujQUkhij^ (20) 



shows that all of the terms involving A cancel out, as they must, leaving the condition 



Iq + riiJi = -uj n l n j h l:j (0). (21) 



3. Photon is emitted with frequency coq 



The photon is emitted at the source with the unperturbed frequency ujq — —g^k^u 1 ^. The 4-velocity of the 
source is u£ = (1, 0, 0, 0) as it has constant spatial coordinate position, so the constraint becomes —g^k* 1 = ujq. 
The source emits the photon at A = X s , so from Eq. (|T5|) this yields 



h = ^nVZ„-(A s ) (22) 

4. Perturbed photon path must hit source worldline somewhere 

The constraint that the perturbed photon trajectory must hit the source worldline somewhere can be written 
as 

X s (X s )=x>= 4 (A.) + 4 (A. ) (23) 

for some A s . To zeroth order we have A s = A s , but there will be a first order correction. Inserting the expression 
(fT6|) for the perturbation of the geodisic gives 

x{ = -LJ X s n j - ^tfSaiXs) + X s Jg . (24) 



Projecting this equation perpendicular to n gives a formula for the perpendicular component of Jq: 



4± = ^{5 ik -rjn k )riS jk (\ s ). (25) 

Here on the right hand side we have replaced X s with A s , which is valid to linear order. Adding to this our 
earlier result for the component of Jq parallel to n in Eqs. (|2~T|) and (|22|) gives 

4 = ^-n j S jk (X s ) (<f fc - n l n k ) ~ ^-n l n^n k 2 jk (X s ) + ^ nV'n fe /i, fe (0). (26) 

D. Perturbation to Observed Frequency 

We calculate the observed photon frequency uj b s — —g fJ ,uk^u^ bs , where u v ohs — (1,0,0,0), and check our result 
against standard formulae for the frequency shift, used in pulsar timing searches for GWs [39(. The observed frequency 
is, from Eqs. {IBJ and (J22J, 

uiobs = fc°(0) = w + Jo - + L ^-n i n j Iij(X s ). (27) 
Using the definition (|15l) . the perturbed redshift is therefore 

/•A s 

z= = — —nn J dXhijfl(X). (28) 

wo 2 J 

For a plane wave traveling in the direction of the unit vector p, we have 

hj = hij(t - p x) = hij [u) X(l + p • n)] , 

giving 



h = 9 h 1 d h 

where 7 = p • n. This gives for the redshift 



1 



Z = ~ 2(1 + 7) " V ^ l ^ A ^ ~~ ' 
which agrees with [3^ | up to a sign, which is an error in their calculation [40l |. 



(29) 



E. Local Proper Reference Frame of Observer 

We must also account for the changes induced in the basis vectors of the observer's local proper reference frame 
due to the presence of the GW. We introduce a set of orthonormal basis vectors e& which are parallel transported 
along the observer's worldline, with eg = u. The parallel transport equation for the spatial vectors gives 



u a e? =u° 



= 0. (30) 



We separate the basis vectors into two pieces, e\ — 8\ + de\, where we assume that the unperturbed basis vectors are 
aligned with the coordinate basis directions. 



Using u — dt, and the connection coefficients (HJ) of the metric , Eq. (|30j) gives us an expression for the perturbation 
to the basis tetrad: 



8e\(t) = --h i3 (t,0) + u,$, 

where u l ~ is a matrix of constants. Now, we observe that e- is an orthonormal set of three- vectors, which gives us six 

j J 

constraints on the constants oji : 

3 

(W + h mn ) (ST + Sef) (d% + Sefj = S ]k . 

This is identically correct to zeroth order; to first order we get Se^ + 5e k ~- + h jk = 0, or, inserting our equation for 
5e, and assuming that hij = hji, we find = —u>ji, i.e. that the constants are antisymmetric in their indices. 
These constants parameterize an arbitrary infinitesimal time-independent rotation. Evaluating now at the detection 
event gives 



Se\ = -\h i3 (0)+ Ui3 . (31) 

For the remainder of this paper we will set to zero the term u> i3 , since it corresponds to a time-independent, unob- 
servable angular deflection. The deflections caused by GWs will be observable because of their time dependence. 

F. Observed Angular Deflection 

We can express the four-momentum of the incoming photon in the above reference frame as 

k a {0) = io obs u a - u obs nhf, (32) 

where S'-piPn = 1, u a is the observer's 4-velocity, and uj b s is given by equation l27l Note that we evaluate all 
quantities at the detection event t = to, x = 0. Plugging in our results for the perturbed 4-momentum and the 
observed frequency, we obtain an equation for the observed direction to the source n° 

fc*(0) = - u n l + ^-n j S jk (X s ) (S ik - n l n k ) - ^nV 'n k T jk {X s ) + \u Q n l n^ n k h 3k (Q) 



= - ( cjo + -£n k n%i(\ s ) ) n> ( 5) - -h)(0) ) . (33) 



We decompose the direction to the source into zeroth and first order pieces as n? = n 3 + Sn^ . The zeroth order 

terms in Eq. ()33j) give us n 3 = n 3 . Plugging this into the first order terms and simplifying, we find the perturbation 
to the source direction 



Sn 1 = - \ n j hij(Q) - -V S jk (X s ) (S lk - n l n k ) - n l n 3 n k h 3k {0) 
Inserting our definition of Sj k , we obtain the solution to the source direction perturbation in Minkowski space 



Stf =-\ n j hij(Q) - 72 l n 3 n k h jk (0) - ^ (S lk - n l n k ) n j 



A., 



A' 



A' 



dX' dX"h jk , (X") + ri dX' dX" (h jk ,i(X") + h kld {X") - hji^X")) 



(34) 



As a check of the calculation, we see that 8n % is orthogonal to n 4 , so that n l + 8n l is a unit vector, as expected. 
We now specialize to the case of a plane wave propagating in the direction of the unit vector p, 

hij(t,x) = hij(t - p • x). 

Using the identity (j2"0)l we can reduce the double integrals in Eq. to single integrals, obtaining 



^' fc(A) -2(l + p.„)^ (A) 



(35) 



Evaluating this explicitly for the plane wave 



hij(t, x) = Re 



gives 



Stf =Re 



i(2 + p- n) 



1-e" 



iQcjQ (1+p-n) A £ 



woA s ri(l + p • n) 

^ _ g— ir2wo(l+p-n)A s 



w A s ft(l + p n) 



<?n k U jk e~ iato 
2(l + p-n) 



1 



2 cjoA s f2(l + p • n) 



^ _ g-jf2tJo(l+p-n)A a 



(36) 



If we define the observed angles (9,<fi) by n* = (sin^ cos <fi, sin sin 0, cos 9), then the observed angular deflections 
are 



e\5n l 



sin 



e ' 



(37) 



where e\ = (cos 9 cos 0, cos 9 sin 0, — sin 9) and e z - = (— sin <j>, cos 0, 0). 

As another check of our calculation, we now compare our result with the coordinate (gauge-dependent) angular 
deflection computed by Yoo et al. [4l|. Starting from our Eq. we disregard the first term, which arises from 

the change from the coordinate basis to the parallel transported orthonormal basis. The remaining terms in Eq. 

give the coordinate angular deflection Sri 1 . Simplifying using the identity (|2T))) and the identity dx' J* dx" f{x") — 
f Q dx'(x ~ x')f(x') gives 



6n l = -inWn*M0) + - £' dX^^^-^n k + ^ (^^) d d (nV/i w )j . (38) 

When combined with Eqs. (|37p . this agrees with Eqs. (13) and (14) of [4l|, specialized to only tensor perturbations, 
up to an overall sign. The sign flip is due to the fact that Ref. [41| uses a convention for the sign of angular deflection, 
explained after their Eq. (16), which is opposite to ours. 



G. The Distant Source Limit 



We now specialize to the limit where the distance uio | A s | to the source is large compared to the wavelength 
~ c f2 _1 of the GWs. As discussed in the Introduction, astrometry is potentially sensitive to waves with a broad 
range of frequencies, extending from the inverse of the observation time (a few years) down to the Hubble frequency. 



Therefore this assumption is a nontrivial limitation on the domain of validity of our analysis. However, for sources at 
cosmological distances (the most interesting case), this assumption is not a significant limitation. 

In this limit, we can neglect the second term in each of the three small square brackets in Eq. (|36|) . giving 



Sn\t,n) = Re 



-itlt 



2(l + p-n) 



-ittt 



(39) 



where we have written t for to- This result agrees with and generalizes a calculation of Pyne et al. [33T ]. We note 
that this same approximation is used in pulsar timing searches for GWs [39j . In that context the approximation is 
essentially always valid, since pulsar distances are large compared to a few light years, and the properties of pulsar 
frequency noise imply that that pulsar timing is only sensitive to GWs with periods of order the observation time, 
and not much lower frequencies, unlike the case for astrometry. 

III. GENERALIZATION TO COSMOLOGICAL SPACETIMES 

Of course, we do not live in Minkowski space. The apparent homogeneity and isotropy of the universe imply that 
our universe has an FRW geometry, with line element: 



ds 2 = g a /3dx a dx l3 = a(r) 2 {-dr 2 + [8 l3 + /iy(r,x)] dx l dx j ) 



(40) 



where t is conformal time, and we specialize to the transverse traceless gauge in which S^hij — S tJ dihjk = 0. To 
translate our calculation in Minkowski spacetime to this new metric, we define an unphysical, conformally related 
metric g a p = a{T)- 2 g a/3 given by 



g aj 3dx a dx^ = —dr 2 + [5ij + ft.y(T, x)] dx l dx-* , 
which has an associated unphysical derivative operator V Q . 



(41) 



A. Stationary Observers are Freely Falling 

As before, it is straightforward to check that observers who are stationary in the coordinates (|4Tfl) are freely falling. 
Therefore we assume as before that the observer and source are stationary: 

xi bs (t)=o, xi(t)=xi 

B. Null Geodesic in the Conformal Metric 

Let us consider a photon traveling from a distant source to us, which follows a null geodesic in the physical metric 
g a /3. Its path is also a null geodesic of the conformally related metric g a p, though it it is not affinely parameterized 
in this metric [42|. Specifically, the physical 4- momentum of the photon must satisfy the geodesic equation 
fc M V \j,k u — 0. If we define a conformally related, unphysical 4-momentum — k^, whose contravariant components 
are then related to those of the physical 4-momentum by 



k" = g^K = a(T) 2 g" u ku = a(r) 2 k v = a(r) 2 fc^, 



(42) 



then we find that 



k^V^K = a{T) 2 k^V^k v . 



(43) 



From |42| we know that for any vector v a , and conformally related derivatives V Q and V Q , we have = Vq,?;^ 

C*2 7 w 7 , where C^ 7 = 2^ Q V^) In a — g Q ^g 7<5 V 's In a. Thus, we find 



k^^h = a{T) 2 k^V^k v + a{r) 2 k^k p fosfr?^ In a - g^g pa V ff In < 

= a{T) 2 k^V f.K + a(r) 2 {k p k p V u In a + k p k v V ^ In a - k v k a V a In a) 

= a(r) 2 fc"V,t, (44) 

where to get the last line we have used that the geodesic is null. Therefore, if fc M is a null geodesic of the physical 
metric g^, then k^ is a null geodesic of the conformally related metric g piy . If A is an affine parameter of the geodesic, 
it will not be an affine parameter of the geodesic in the unphysical metric. The affine parameter A in the unphysical 
metric is related to A by 

dX a(r(A)) 2 ' [ ' 

C. Parallel Transport of Basis Vectors in FRW Background Spacetime 

We next investigate the parallel transport of the observer's basis tetrad in a FRW background spacetime. From 
the form (|40p of the metric, we anticipate that the basis vectors must scale as a -1 to remain normalized. Thus, we 
will define the basis vectors and their perturbations as 



e-. 
j 



+ (46) 



The relevant connection coefficients are 



n k = -Jl+ l -8™h mkfi . (47) 



The parallel transport equation (|30|) for the spatial basis vectors gives us 



a,ei+r* fc ej = 0. (48) 



Plugging in the connection coefficients (|4T)l and the basis vector expansion (|4"6| . we get 



doSel + \h m h m] = 0, (49) 



the same equation as before. The solution, as before, will be 



SeUt) = --h)(t). (50) 



D. Generalization of Angular Deflection Computation 

We parametrize the photon path in the background spacetime by 

t(C)=to+ Wo C, x\0 = -CLo a n\ (51) 

where Q is an affine parameter of the unphysical metric (|41j) (denoted A above). From the decomposition (|32[) . the 
observed source direction is 



We rewrite all the quantities in this expression in terms of their conformally transformed versions 



e~ = ae~. , 

3 3 ' 



9 a p = a 2 g a p, k a =a 2 k a , u a = an 
which are the quantities that are used in the Minkowski spacetime calculation of Sec. [TXJ This gives 



(53) 



(54) 



the same expression as in Minkowski spacetime. Therefore, the final result is the same expression ([Ml) as before, 
except that it is written in terms of the non-afhne parameter £: 



Sri 1 =-\ n J %(0) - n i n j n k h jk {0) - y- (S ik - n i n k ) ti? 

2 l°d£ t dC'h jk , Q (C') + n l ( S d(' [ C dC'ihjk^O + hkijiO-hjiMC")) 



n 



o 



(i 







(55) 



E. 



The Distant Source Limit 



We now specialize again to the limit where the distance to the source is large compared to the wavelength ~ cQ^ 1 
of the GWs. We also assume that the wavelength c f2 _1 is small compared to the horizon scale, but we allow the 
sources to be at cosmological distances. 

Starting from Eq. (|55|) and paralleling the derivation of Eq. (|3"5j) we obtain 



6n z (T , n) = ^SikUjhjkiO) + / d( [h jk (Q - h jk {0)} + 



d( ( ^^ - Tijnihji,k(0 ) , ( r > ( >! 



where Si 



UiUk- Now the wave equation satisfied by the metric perturbation is 



d 2 + 2-^d T - V 2 J hijfax) = 0, 



and plane wave solutions are of the form 



Mr.x) =Re{U ij e inpx qn(T)} 



where the mode function g$i satisfies 



1- 
a 



2-q' n + Q 2 qn = 0. 



(57) 



We now evaluate the angular deflection (|56l) for such a plane wave, in the limit where e = a' /(Via) <C 1, i.e. the 
limit where the wavelength ~ a/ft of the GW is much smaller than the the horizon scale ~ a 2 /a' . In the second term 
in (1561) . the term hj k (C) is rapidly oscillating, and so its integral can be neglected compared to the integral of hjk(0)] 
corrections will be suppressed by powers of e. In the third term in the integrand is rapidly oscillating, and so 
the integral will be dominated by contributions near the endpoints, up to 0(e) corrections. However the integrand 
vanishes at £ = £ s , and thus the integral is dominated by the region near C = 0. In that region we can use the leading 
order WKB approximation to the mode function solution of (|57| , 



and to a good approximation we can replace a(r) by a(ro). Thus we see that the same answer is obtained for distant 
sources as in our Minkowski spacetime calculation, even for sources at cosmological distances. From Eq. (|39|) we 
obtain 



Sn l (T , n) = 2(i + p. n ) ^ fc (°)"j nfc - ^MOK ^ 



for plane waves in the direction p. 



IV. CALCULATION OF ANGULAR DEFLECTION CORRELATION FUNCTION 



Now that we have calculated the deflection of the observed direction to a distant source due to an arbitrary metric 
perturbation hy, we would like to determine the properties of the deflection produced by a SGWB, such as that 
produced by inflation. 



A. Description of SGWB as a Random Process 



In the distant source limit, the angular deflection (|58| depends only on the GW field hij evaluated at the location 
of the observer for each direction of propagation p. Moreover, we have restricted attention to modes with wavelengths 
short compared to the Hubble time. Therefore, it is sufficient to use a fiat spacetime mode expansion to describe the 
stochastic background. This expansion is (see, e.g. Ref. [27| ) 



hy(x,t) = £ r $ I d " n p W/) e 2 ^-*-^ + c.c, (59) 

A=+,x Jo J 

where / and p are the frequency and direction of propagation of individual GW modes, h Ap are the stochastic 
amplitudes of modes with polarization A and direction p, and the polarization tensors efj' p are normalized such that 

We will assume that hij(x, t) is a Gaussian random process, as it is likely to be the sum of a large number of random 
processes. We also assume that it is zero-mean and stationary. It follows that the mode amplitudes h,A P {f) satisfy 



(h Ap (f)h Bp '(f')) -0, 

(h Ap (f) h Bpl {fY) = ^T^P *U ~ /') Sab <5 2 (P, p') (60) 



for /, /' > 0, where Hq is the Hubble parameter and <5 2 (p, p') is the delta function on the unit sphere (see, e.g., j27j]'). 

Since the angular deflection 5n(n, t) depends linearly on the metric perturbation, it will also be a stationary, zero- 
mean, Gaussian random process, whose statistical properties are determined by its two point correlation function 
(SrtSni). Specializing our expression (155)) for the angular deflection to the form (f59"|) of the metric perturbation, we 
find 

Sn l (n,t)= V / df ( d 2 tt p h Ap {f)e- 2mit K lkl {n,v)e^ + c.c., (61) 
^_ „ Jo J 



where 



-o , \ 1 f[ni+Pi]n k ni , 

K ik i(n,p) = - —— n k S u . (62) 



2 V 1+p n 



B. Power Spectrum of the Astrometric Deflections of the SGWB 

So, we need only evaluate the two-point correlation function to gain full knowledge of the statistical properties of 
the angular deflection due to the SGWB. Writing out this quantity explicitly using Eq. ([M]) . 



{8n\n,t)8ni{n',t')) = V / dfdf f d 2 fi p d 2 O p V \h Ap (f) n ikl (n,p) e% p + c.c. 

A,b=+,x Jo J \ l J 

x [W(/T ^ lft ' ^> s (n',p') (ef/)* + c.c] V (63) 



The average, which is an average over ensembles, acts only on the stochastic amplitudes h,A P - Using the mode 2 point 
function (|60|) in Eq. (|63| . we get the simplified result 

/•CO o Til 

(5n\n,t)5n>(n',t'))= / df ^/-^(/Je"^'-^-^, n') + c.c, (64) 



(i 



where we have defined 



H l3 (n,n')= [ d 2 n p 1l ikl (n,p)e% p K jrs (n',p) (e^ p )*. (65) 

4_ i "\/ 



A= + , 



C. Basis Tensors and their symmetries 

We simplify the expression (|65[) for Hij further using the identity 

E 4' p (4r)*=iPw> (66) 

A=+,x 

where Piju is the projection tensor onto the space of traceless symmetric tensors orthogonal to p, given by 

2-Pyfcz = S ik Sji + SuSjk - SijSki +PiPjPkPi - SikPjPi - SjiPiPk - SuPjPk - SjkPiPi + SijPkPi + hiPiPj- (67) 
This gives 

H ij (n,n')=2 J d 2 n p n i ki(n,p)P k i rs n jrs (n' ,p). (68) 

Noting that the correlation function (|64j) is perpendicular to n on its first index and n' on its second, we can 
decompose it onto a basis of tensors with this property: 

Hij{n, n') = a(n, n'^A/ + /?(n, n!)A l C J + 7 (n, n'^A, + cr(n, n!)B t C v (69) 
for some scalar functions a, /3, 7 and <r. Here we have defined 

A = n x n', B = n x A, C = -n' x A. (70) 

We can deduce from Eq. that iJy (n,n')* = H 3l (n',n). Noting that Aj(n',n) = — Af(n,n'), and 5j(n',n) = 

— Cj(n, n'), this symmetry applied to the expansion (iMl) gives 



a(n, n')* = a(n',n), cx(n, n')* = o-(n', n), /3(n, n')* = 7(n', n). 



We see from Eq. (|65j) that ffy transforms as tensor under rotations. This implies that the functions a, (3, 7 and a must 
be invariant under rotations, and can only depend on the angle 9 between n and n'. Thus, a(n, n') = a(n', n) = a(8) 
and so forth, so a and a must be real. 

Next, we note that the expression for ffy(n, n') is invariant under the parity transformation n — > — n and 
n' — > — n'. Looking then at the basis tensors, we see that A is invariant under this transformation, while B and C 
change sign. Thus, in order to insure that ffy is invariant, it can only have terms multiplying AiAj and BiCj, so 
0(0) =0 = 7(8). 

Having taken the symmetries of the problem into consideration, we have found fiy to be of the form 



#y(n,n') = a{Q)A i A j + a{Q)B i C j 



(71) 



D. Solving the General Integral 

We can evaluate the coefficients in the expansion (|7Tjl of by contracting it with the basis tensors: 

.l'.W/,, = sin 4 (9)a(9), ll'C'll,, = sin 4 (9)c7(9). 
Rewriting these using Eq. (|68j). we find 



a(Q) = , * [ d 2 n p A t n m (n,p)P klrs A^TZ jrs (n',p) 
sm (9 J 



(72) 



a(G) 



sin 4 (9) 



d 2 n p B l K m (n, p)P klrs C j 1l jrs (n', p) 



(73) 



To simplify the calculation, we define the quantities k = n • p, k' = n' • p, A = n • n', fj, = A • p, which satisfy 
/i 2 + A 2 + k 2 + k' 2 = 1 + 2Xkk'. Using these definitions and the definition (p)2"j) of 7^^, we can write 



A i n m {^p) = \n k (^- 

i \ 1 + K 



A) 



1 / k' + A A 
B 7^ fci (n,p) = -n fc I - - - - m + n, I , 

We can then rewrite our expressions for a and a 



1 



A J 7^ irs (n',p) = -n 



1 + k' 



A, 



C j K jrs (n',p) 



1 , / k + A , 



1 + K 



a ^ = J-^T / d2 ^ 2P ^ s n k A t ) < (-^L - A s ) , 

4sm (9) J \1 + k / + K / 



(74) 



a(9) 



4 sin 4 (9) 



, << > > ,-, / k' + A A . / k + A , 
<i ll p 2P klrs n k [ — 1 , + ) n r | -— — 7 n s + n. 



1 + K 



1 + k' 



(75) 



Let's define two new variables z^ 2 = (1 — k 2 ), v 12 = (1 — k' 2 ). Applying the definition (|57| of the projection tensor 
Pkirsi we can calculate the necessary contractions of P k i rs for a: 



2P Wrs n fe Ai?<A = (A - kk') (1 - A 2 - 
2P klrs n k Ain' r n' s = \i (kk' 2 - 2Xk' + k) , 



2P klrs n k niri r A s = (k'k 2 - 2Ak + k') , 
2P k i rs n k nin' r n' s = v 2 v' 2 - 2/i 2 , 



(76) 



and for tr: 



2P k i rs n k n' l n' r n s 
2P k [ rs n k nin' r n' s 



2P klrs n k n l n' r n s 
2P k [ rs n k n' l n' r n' s 



= v 2 (A — kk') , 
= v 12 (A - kk') . 



(77) 



Plugging these back into Eqs. ([74]) and (1751) and simplifying, we find 



(A- KK )(1-A)-M (1 + A) + (i + k)( i + k0 



-<r(0). 



Noticing that we can do the integrals J d 2 fl p fi 2 = 4^ sin 2 8 and J d 2 fl p KK' — 4^ cos 8, but that the last term is more 
complicated, we find 

a(e ) = - g (e) = ; (c0Se - 1) + -L 7 , /^O p ^ (A + l (A + / f ) . (78) 



3 sin 2 8 2 sin 4 8 7 P (1 + /e)(l + k') 

We can reduce the two dimensional integral (178[) to a one dimensional integral by parameterizing p in spherical 
polar coordinates 9 p and <p p , choosing n = (0, sin(6/2), cos(6/2)) and n' = (0, — sin(8/2), cos(6/2)) and integrating 
over <pp. This gives 

a(8) = -(7(0) = - v . 2 ; , + , 5 / d0 p sm6 p {sin 2 P + 8cos(8/2) [cos6> p + cos(8/2)] [g(0 p , 6) - 1]} , 



3 sin 2 (8) 2 sin 2 8 Jo 

(79) 

where 

|cosfl p + cos(8/2)| 

g( ^' 0) - [f + cos^cos(8/2)]- (80) 
We perform the integral over 9 p , and find the final form of the function a(8) 

7T 327T 

a(0) = -a(0) = —3- (7 cos 8 - 5) - — In (sin(8/2)) sin 6 (8/2). (81) 
3sm 8 sin 8 

A plot of the function a(8) is shown in Fig. [5J 

To summarize, we have now completed the calculation of the angular deflection correlation function. The final 
answer is given by Eq. (|64|) . with ffy (n, n') given from Eqs. (fTTj) and (j8Tj) as 

tfy (n, n') = a(0) (A A, - B L C 3 ) . (82) 
Here the vectors A, B and C are defined by Eqs. (|70|) . and a(8) is given by Eq. (j8Tj) . 

E. Special Case: Coincidence 

As a check of our calculation, we can solve for the two-point correlation function exactly in the case that n = n'. 
Using Eqs. (1B2"|) , (p7|) and (frJSj) . the integral simplifies to 

Hij(n, n) = ~ y d 2 O p [l - (p • n) 2 )] (% - n.nj) . 
We can solve this integral analytically, getting 

Hij(n, n) = — (<5# - mnj) . (83) 



This corresponds to the limit of a(Q)(A i A j - &&) as n ^ n', with a(0) = 2tt/(30 2 ) + ©(e^ 1 ) from Eq. (JgTJ. 
Inserting the coincidence limit (|53")) into the correlation function yields the formula © for the total rms angular 
fluctuations discussed in the introduction. 
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FIG. 2: Here we plot the function a(G), the coefficient of Hij(n, n') as shown in Eq. {82) , as a function of the angle O between 
n and n'. 



V. SPECTRUM OF ANGULAR DEFLECTION FLUCTUATIONS 

A. Overview 

In the previous section we computed the correlation function (<5n J (n, tySni (n', t')) as a function of the unit vectors n 
and n'. However for many purposes it is more useful to perform a multipole decomposition of the angular deflection, 
and to compute the spectrum of fluctuations on different angular scales I , as is done with cosmic microwave background 
anisotropies. We decompose 5n(n, t) as 

<5n(n,t) = J2 Sn Elm(t)Yf m (n) + «n Blm (t) Y,«(n), (84) 

lm 

where Y^ and Y^ are the electric- and magnetic-type transverse vector spherical harmonics defined by 

Y? m (n) = (1(1 + l))- 1 /2 V y (m (n), Yf m (n) = (1(1 + l))-^ 2 (n x V)Y lm (n). (85) 
We will show in this section that the statistical properties of the coefficients are given by 

(8n Qlm {t)6n Qn - m '(t'T)=5 QQI 8 lv 5 mm , / df cos [2tt/ (t - t')]S Ql (f) (86) 

Jo 

for Q, Q' = E or B, for some spectrum Sgi(f), a function of frequency / and of angular scale I. The formula 
(|86| shows that different multipoles of the angular deflection are statistically independent, as required by spherical 
symmetry of the stochastic background. Also the electric-type and magnetic-type fluctuations are uncorrelated, as 
required by parity invariance of the stochastic background (see below). 



The spectrum Sqi(J) is given by 



%(/) = 



An 
2Z + 1 



9 2 

rm.s 



/ 



(87) 



Here # 2 ms is the total rms angular fluctuation squared, given by Eq. (j3|) in the introduction. The function <r(f) 
describes how the power is distributed in frequency. It is the same for all multipoles, is normalized according to 
J d( In /) cr(f) = 1, and is given explicitly by Eq. (|6]) in the introduction. The quantities and gs are the fraction of 
the total power in electric-type and magnetic-type fluctuations, and are gE = gB = 1/2. Finally the angular spectra 
af E and af B describe the dependence on angular scale, which is the same for all frequencies. They are normalized 
according to 



E 

i=i 



*? Q = i, 



and are the same for E and B modes, af E = af B . This spectrum is plotted in Fig Q] and the first 10 values are listed 
in Table HI We note that these coefficients are well fit by the power law af E — 32.34 Z~ 4 - 921 . 

Before proceeding with the derivation of the spectrum (l86l) . we first derive from (|86|) the expression ((4]) discussed 
in the introduction for the total fluctuation power. Squaring the expansion taking an expected value, and then 
using ([55]) gives 



Qlrn Q'l'm' 

I 



Ql J0 J m=-l 



9 2 
rms 



Ait 
21 + 1 



9Q a i 



QQ 



Using Unsold's theorem for vector spherical harmonics, 



(89) 



21 + 1 
4tt ' 



gives 



/•OO 

Ql Jo 



2 

rms 



9q «r 



(90) 



which reduces to Eq. Q. Note that using the normalization conventions for and a(f) now gives (<5n(n,t) 2 ) 
@rms(9E + gB) = Orms' showing consistency of the definitions. 



B. 



Derivation 



We now turn to a derivation of the spectrum (|87j) . First we note that the vector spherical harmonics are transverse 
in the sense that Y^(n) • n = for Q — E, B, and are orthogonal in the sense that 



d*n n Y l ( i t (n)Y l ?J(n)=6 QQ ,6 w 6 n 



Using this orthogonality property, we can extract the coefficients of the expansion (|84 



5n 



Qlrn 



(t)= J d 2 fi n ^(n,^"(n). 



Thus we can write for the correlation function between two of these coefficients 



{5n Qlm {t)5n Q , Vm ,{t>)*) = J d a n n d 2 n n ,y«*(n)^, i (n , )<«n i (n,i)«n J (n / ,t')) J (91) 
or more explicitly, using Eq. 

(6n Qlm (t)5n Qn , m ,(t')*) = J dfcos[2nf(t - if)] -^U. C QlmQ , Vm ,, (92) 

where 

C QlmQ n' m > = j d 2 0„d 2 ^y«*(n)y^(n')^(n,n'). (93) 

We now argue that the EB cross-correlation vanishes. From Eq. we see that YjL(n) has the same parity under 
n — > n as Yj m (n), while the parity of Y B n (n) is opposite. From section HV CI above Hij(n,n') is invariant under 
both n — > — n and n' — > — n'. Thus, if Q = E, Q' = B in Eq. (j9"2"]l . the integral will be symmetric under n — > — n but 
antisymmetric under n' — > — n', causing the integral over <i 2 f2 n < to vanish. Therefore, EB cross correlations vanish, 
and we need only calculate the EE and BB correlation functions. 



1. EE correlation 



Inserting the definition (1551 of the electric vector spherical harmonics and the formula (|52")l for Hy into Eq. (|9"3")l 
and integrating by parts, we obtain 

C E lmE>Vm' = J d^d 2 ^, Y* m {n)Y v m , (n')/3 BB (9) , (94) 
where the function j3 EE is given by 



f3 EE (Q) = V,V;- [Hy (n, n')] = V<V$ {a(0) [A^- - B^-]} . (95) 

Here Vi and V^- denote normal three dimensional derivatives with respect to x and x', where n = x/|x| and n' = 
x'/|x'|. Integration by parts on the unit sphere of this derivative operator is valid as long as the radial component of 
the integrand vanishes, from the identity ViU* — d r v r + 2v r /r + S7av a , where V^i denotes a covariant derivative on 
the unit sphere. It can be checked that the radial components do vanish in the above computation. 

Next, we expand the function f3 EE in terms of Legendre polynomials, and use the spherical harmonic addition 
theorem, which gives 

i 

= E ^C^W^tnT (96) 

Ira 

Inserting this into Eq. (|94p and using the orthogonality of spherical harmonics gives 



CElmE'l'm' — <>U'0~rnm' Jq ^ ~ ^ Pf E ■ (97) 

Inserting this into Eq. (p?2"j) now yields the correlation function given by Eqs. (jH6"| and (1571) . and using the definitions 
(U) and ^ of &l ms and a(f) allows us to read off the electric multipole spectrum 



(98) 



We will show below that gE — 1/2. 

It remains to explicitly evaluate the function (3 EE (Q) defined in Eq. (|95|) and evaluate its expansion coefficients. 
We have 

(3 EE {G) = V 4 V; [a{Q)T ij ] = [V,:V;.a(8)] T l] + [V<a(0)] (V^) + [V<a(0)] (V<T«) + a(0) (ViV^) , (99) 

where we have defined T ij = (AM^n, n') — B i C j (n,n')) . Using A 4 = e ijk njn' k , B l = (n • rr> 1 - n'\ C* 1 = 
(n • n')n n — n 4 , we can write the tensor T ZJ in Cartesian coordinates as 

T ij = e ^ e ^ nkn ' inrn ' s _ (( n • nV - n H ) ((n • n')n' j - n j ) . 
Using ViTij = Sij — riifij, V^n^ = &y — n^n^, Vjn^ = V^n' 3 = 0, and V;e lJ = Vje yfe = 0, we calculate the derivatives 

ViT y = (1 - 3(n • n')) ((n ■ n>' J ' - re*) , V'jT** = (1 - 3(n ■ n')) ((n • n')n* - n H ) , 

ViVjT^ = - 9(n • n') 2 + 2(n • n') + 3. (100) 

For the gradients of a, we use the fact that cos(O) = n • n', so that — sin(0)Vi0 = v! i — (n • n')nj, and similarly 
for V'y Thus, we find 



sm(9) J sm(fc>) 

Sij - riiUj - n^n'j + (n • n')mn^ cos(9) [n- - (n • n')m] [rij - (n • n')n'j 



ViV' i a(Q) = a'(e) 



■sin(0) -sin 3 (0) 

,„ /Q ^ K - ( n ' n >»] ["j - ( n ' n 'K] 
sin 2 (6) 

Plugging Eqs. (1TUU|) and ([TUB into Eq. JSUJ, we get 



a"(6) — — 21L — • ( 101 ) 



{3 EE {Q) = [-9 cos 2 (6) + 2 cos(6) + 3] a(0) - sin 2 (6)a"(e) 

+ [l-6cos(0)]sin(0)a'(e). (102) 

Next, we insert the expression (j8Tj) for a(0) to obtain 

^ E (6) = ^4 + (1 - cos 9) {12 In [sin(0/2)] - 1} ) . (103) 

We numerically compute the coefficients f3 EE of the Legendre polynomial expansion (|UU[) of f3 EE (0), and from them 
compute af E using Eq. (jUgj) . The result is plotted in Fig. [JJ and tabulated in table UJ 

2. BB correlation 

We now calculate the BB correlation in a similar manner to the EE case above. Inserting into Eq. (|93|) the 
definition ([55)1 of magnetic vector spherical harmonics and integrating by parts, we find 

CBimBi'm* = jjj^-fi J d 2 n n d 2 Q n , Y? m (n)Y Vm . {n')p BB (Q) , 



where 



(3 BB (<d) = V ; V; [e ik ie jmp n k n' m a{Q)T i3 ] . (104) 
As before, we can derive from here the form (|86|) and (|87|) of the spectrum, with af B given by 

„ „BB ^ qBB 

9Bai = i^(7TT) A ■ 

We now show that (3 BB (0) = f3 EE (<d), from which it follows that g E = 9b = 1/2 and that af E = af B . To see 
this we evaluate the cross products in (|104[) using nxA = B,nxB = —A, n' x C = A. This gives 

and using the definitions ^ and (HH1 of /3 EE and (3 BB , it follows that (3 BB = (3 EE . 
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Appendix A: Non Local Dependence of Deflection Angle on Gravitational Waves 

It is sometimes claimed in the literature that the deflection angle, in the distance source limit, depends only on 
the GWs in the vicinity of the source and the observer. A similar claim is often made for the frequency perturbation 
caused by GWs which is targeted in pulsar timing searches. Strictly speaking, these claims are not true: it is possible 
to have a nonzero deflection when the GW field vanishes in a neighborhood of the source and of the observer, even 
in the distant source limit. However, this type of circumstance requires a considerable fine tuning, so the claims are 
colloquially valid. 

To see this, it is sufficient to consider the simple model of a scalar field h(t, x) in 1 + 1 dimensions, obeying 
(<9j — h = 0. A functional of h(t,x) that is qualitatively similar to the deflection angle formula (|64[) is 

A9 = c R h R (t - .T obs ) + c L h L (t + x ohs ). (Al) 
where we have decomposed the field into left-moving and right-moving pieces, 



h(t,x) = h R (t-x)+h L (t + x), (A2) 

and c R and cl are fixed coefficients with cl 7^ c R . In each sector (right moving and left moving), the quantity (|A1[) 
depends only on the field evaluated at the location of the observer x = x bs- However, the sum does not. If we specify 
the field in terms of its initial data h(t, x) and h(t, x) at time t, t, and specialize to initial data of compact support, 
we get 

1 1 r x "»s 

A9(t) = -(c L +c R )h(t,x obs ) + -(c L -c R ) dxh(t,x). (A3) 
" " J -00 

It is clearly possible to choose h and h to vanish in a neighborhood of x — x i, s and still have A9 7^ 0, for cl 7^ c R . 



